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ABSTEACT 

The  problem  ia  treated  of  an  elllptlcally  polarized  plane  electro- 
magnetic wave  incident  obliquely  on  a  cont inuoiiBly  stratified  anisotropic 
scattering  mediiam.  The  medium  occupies  a  region  of  space  "between  two  paral- 
lel plBjies  (O  ^  e  ^  a  say)  and  is  characterized  by  a  Tnagnetic  suscoptluilluy 
which  is  zero  and  an  electric  susceptibility  which  is  a  continuous  matrix 
function  of  z.  It  is  required  to  obtain  the  reflected  and  transmitted  plane 
waves  on  the  incident  and  far  sides  of  the  medium,  respectively. 

The  problem  is  attacked  by  means  of  the  Schwinger  integral  equation- 
variational  technique.  An  integral  equation  is  derived  for  the  electric  field 
vector  within  the  scattering  medium.  Next,  the  problem  is  reduced  to  comput- 
ing the  elements  of  a  certain  •reflection  matrix ■  and  •transmission  matrix", 
A  set  of  variational  expressions  for  these  matrix  elements  is  then  constructed 
by  msfms  of  simple  considerations  in  the  theory  of  abstract  linear  spaces. 
The  variational  expressions  have  the  form  of  functionals  depending  on  two 
independent  vector  functions;  their  Euler  equations  are   the  given  integral 
equation  and  a  certain  %uiJoint  ■  integral  equation,  finally,  a  reciprocity 
theorem  is  proved,  based  on  a  symmetry  property  of  the  variational  es^ressions. 
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1.  Introduction 


The  use  of  variational  principles  for  the  formulation  of  prob- 
lems in  mathematical  physics  is,  of  course,  verj'  old;  one  thinks  of 
Haoilton's  principle,  the  principle  of  least  action,  and  the  various 
kinds  of  mininrum  energy  principles  that  have  been  employed  to  character- 
ize a  multitude  of  physical  situations.  The  attractiveness  of  thase 
principles  is  ascribable  to  a  feature  common  to  all  of  them,  which  is 
that  they  afford  an  extremely  concise  description  of  the  phenomena  under 
investigation.  Moreover,  in  a  great  many  problems  the  variational  prin- 
ciples lead  very  naturally  to  practical  methods  of  estimating  numerically 
the  qijantities  it  is  desired  to  calculate. 

In  recent  years,  the  existing  wide  interest  in  variational  pro- 
cedures has  received  a  powerful  stimulus  from  the  discovery  by  Julian 
Schwinger  that  a  broad  class  of  problems  in  the  diffraction  or  scatter- 
ing  of  radiation  by  obstacles  can  be  formulated  variationally.  These 
problems  are  of  the  type  in  which  a  known  field  of  radiation  impinges  on 
some  kind  of  obstacle  and  it  is  required  to  determine,  at  least  approxi- 
mately, the  spatial  distribution  of  the  intensity  and  phase  of  the  radia- 
tion scattered  by  the  obstacle  at  distances  from  the  latter  which  are 
large  compared  to  a  wavelength;  in  mathematical  terms,  one  desires  the 
asymptotic  form  of  the  scattered  wave  at  infinity.  Schwinger »8  attack 
on  this  problem  consiPts  first  in  deriving  an  integral  equation  for  the 
unknown  field  within  or  on  the  obstacle.  Then  a  set  of  functional  ex- 
preasioas  is  constructed  which  are  stationary  with  respect  to  small  vari- 
ations in  the  function  on  which  they  depend,  in  the  neighborhood  of  solu- 
tions to  the  integral  equation,  ifi,,  the  latter  is  the  Euler  equation  for 
the  functionals.  Tinally,  it  is  shown  that  the  stationary  values  of  these 
f^onctionals  are  precisely  the  values  of  those  numerical  parameters  which 
characterize  the  scattered  wave  at  infinity. 

It  is  this  last  fact  which  is  chiefly  responsible  for  the  prac- 
tical usefulness  of  Schwinger 's  method  since  it  enables  one  to  obtain 
approximations  to  the  far  scattered  field  without  becoming  unduly  involved 
in  the  details  of  the  scattering  process;  one  substitutes  trial  fields 

into  the  variational  functionals,  and  if  these  fields  have  been  skillfully 
*  Hulthe'n,  L.  (K.  Pys.  Sails. Land  rorhand.vol.l^,*21  I9W+)  and  Oinzburg, 
V,L.  (J.Phys.USSB  vol.  7,  p. 289,  19^3)  have  used  alternative  variational 
,v  principles,  but  we  shall  not  deal  with  them  here. 


mmimmmmmmmm^t'^mmrmm^^mmrmmmmmm^m^^mmm^^m 
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chosen  or  perhaps  derived  "bj   iteration  from  the  integral  equation,  then  the 
stationary  character  of  the  fionctioiials  insures  a  determination  of  the 
scattered  wave  parameters  sufficiently  accurate  for  many  applications. 

Schwinger'a  procedure  has  been  applied  very  extensively  "by  him 
to  the  problem  of  the  effect  of  discontinuities  in  electromagnetic  wave- 
guides ;  particular  cases  of  this  problem  have  also  been  examined  by 

2 
Miles  using  Schwinger's  method.  The  question  of  the  diffraction  of 

scalar  waves  by  a  perfectly  reflecting  plane  screen  and  the  complementary 
question  of  diffraction  by  a  plane  obstacle  have  been  dealt  with  by  Levine 
and  Schwinger  •  '  by  means  of  the  variational  formalism,  while  Miles  has 

generalized  this  approach  to  the  aperture  problem  so  as  to  include  electro- 

7 
magnetic  waves.  Sollfrey  has  employed  Schwinger 'a  procedure  to  formulate 

the  problem  of  the  diffraction  of  scalar  waves  by  an  arbitrarily  shaped 

obstacle  and  has  thereby  obtained  a  formula  for  the  scattering  cross-section. 

The  diffraction  effect  of  an  infinite  plane  grating  has  been  treated  by 

8 
Shmoys  through  the  use  of  Schwinger'e  method;  the  grating  consists  of  in- 
finitely many  parallel  wires  of  arbitrary  cross  section,  and  the  waves  are 

9  ' 

electromagnetic.  Pinally,  diTrancia  has  considered  the  case  of  a  plane 

scalar  wave  obliquely  incident  upon  a  plane  stratified  isotropic  medium, 
emd  has  applied  Schwinger '3  procedure  to  obtain  a  variational  expression 
for  the  reflection  coefficient. 

In  the  present  paper,  we  shall  show  how  Schwinger »3  method  can 
be  adapted  to  a  problem  which  is  more  complicated  in  structure  than  any  of 
the  aforementioned  ones.  Essentially,  our  problem  is  a  generalization  of 
the  one  considered  by  diPrancia;  we  shall  be  concerned  with  the  case  of 
an  elliptically  polarized  plane  electromagnetic  wave  with  harmonic  tine 
dependence,  obliquely  incident  upon  a  continuously  stratified,  anisotropic, 
scattering  medium.  The  medium  occupies  a  region  of  space  between  two  paral- 
lel planes  to  which  the  planes  of  stratification  Are  also  parallel.  We  seek 

1.  Saxon,  David  S. ,  Notes  on  Lectures  by  Julian  Schwinger:  Discontinuities 
in  Vfeveguides ,  MIT  Radiation  Lab.  Report,  19^5. 

2.  Miles,  J.  W.,  Proc.  lEE,  vol.  3^,  p.  728,  October,  19^. 

3.  Levine,  H.  and  Schwinger,  J.,  Phys.  Rev.,  vol.  7^,  p.  958,  Oct.  15,  19^. 
k,   Levine,  H.  and  Schwinger,  J.,  Phys.  Rev.  vol.  75,  p.  1^23.1608,  May  1,19^+9. 

5,  Levine,  H, ,  Joum.  Acoust.  Soo.,  vol.  22,  p.  ^,  January,  1950. 

6,  Miles,  J.  W.,  Jotimal  of  Appl.Phys.vol.20,p.  760,  A'ugust,  19^. 

7,  Sollfrey,  W. ,  The  Variational  Solution  of  Scattering  Problems,  Research 
Report  EM-11,  Math. Research  Oroup,  M.Y.U. , Contract  No.AP-19(l22M2 

8,  Shmoys,  J.,  Diffraction  of  Electromagnetic  V/aves  by  a  Plane  Wire  Grating, 
Research.  Report  EI4-18,  Math.  Research  Group, N.Y.U. Contract  No. A?-19 (122)^2. 

9,  di  Prancia,  0.  Toraldo,  II  Nuovo  Cimento,vol.7,p.255iMay  7tl950(ln  Italian). 


-  3- 

the  intensity  and  polarization  of  the  plane  wave  reflected  "back  into  the 
region  from  which  the  incident  wave  arrives,  as  well  as  the  corresponding 
properties  of  the  transmitted  plane  wave  on  the  far  side  of  the  obstacle. 
This  problem  is  of  great  practical  importance  in  connection  with  the  effect 
of  the  earth's  ionospheric  region  on  the  propagrtlon  of  radio  waves  since 
the  planar  model  of  the  ionosphere  can  be  satisfactorily  represented  by 
the  kind  of  inhomogeneous  anisotropic  medivun  Jxist  described.   Indeed,  be- 
cause of  its  significance  for  radio-ionosphere  work,  the  problem  has 
received  some  attention  at  the  hands  of  applied  mathematicians  nsing  more 
conventional  methods  of  analysis  *   *   ,  and  it  is  one  of  the  objects  of 
this  paper  to  demonstrate  that  the  Schwinger  integred  equation-variational 
technique  provides  a  more  succinct  treatment  of  the  problem  than  the  methods 
hitherto  employed  upon  it. 

Another,  and  quite  distinct  point  which  we  shall  make  is  that  it 
proves  highly  effective  to  employ  certain  elementary  considerations  in  the 
theory  of  abstract  linear  spaces,  for  the  purpose  of  assisting  in  the  vari- 
ational formulation.   It  will  be  clearly  seen  how  the  construction  of  the 
variational  functionals  is  expedited  by  the  introduction  of  these  simple 
abstract  considerations. 

The  arrangeaent  of  the  material  in  .this  paper  is  as  follows:  in 
Section  2  we  shall  give  a  precise  matheniatical  statement  of  our  problem. 
In  Section  3,  we  derive  an  exact  vector  integral  equation  which  is  satis- 
fied by  the  unknown  electric  field  within  the  scattering  medium.  Tollow- 
ing  this,  a  certain  ••reflection  matrix"  R  and  •transmission  matrix"  T  will 
be  defined  having  the  property  that  they  operate  on  the  vector  amplitude 
of  the  incident  electric  field  to  yield  the  amplitudes  of  the  reflected 
and  transmitted  electric  field,  respectively;  hence,  our  problem  reduces 
to  computing  the  elements  of  R  and  T.  (If  we  know  the  reflected  and  trans- 
mitted electric  fields,  we  can  obtain  th6  corresponding  magnetic  fields     — 
from  Ifexwell's  equations.)   In  Section  5»  we  shall  state  and  prove  Schwinger's 
fundamental  variational  principle  in  the  fora  of  a  theorem  on  the  station- 
ary property  of  a  certain  ratio  of  inner  products  in  an  abstract  linear 
space.  Each  of  the  elements  of  the  matrices  II  and  T  will  then  be  proved 


10.  Rydbeck,  O.E.H, ,  On  the  Propagation  of  Radio  Waves,  Trans, Chalmers  U.  of 
Tech..  Ho.  dk,   19^. 

11.  WillDes,  M.  y.,Proc.  Royal  Soc.  A  vol.  189,  p.  130,  19^7. 

12.  Bremmer,  H., Philips  Research  Rep.  vol.  4,  p.  189,  19^9. 
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to  eq\»l  the  stationary  value  of  a  corresponding  functional  expression. 
The  form  of  this  erpression  will  "be  deduced  by  setting  up  an  identifica- 
tion between  the  elements  of  the  abstract  space  and  the  vector  f\mction8 
of  the  original  problem  and  inserting  the  latter  into  the  inner  product 
expression.  The  variational  functionals  all  depend  on  two  Independent 
vector  functions,  and  will  be  shown  to  become  stationary  when  the  vector 
functions  are  chosen  respectively  to  be  the  solutions  to  the  integral 
equation  and  to  a  certain  %xLJoint"  Integral  equation. (Although  we  have 
not  hitherto  mentioned  this  adjoint  equation,  its  introduction  is  an 
essential  feature  of  the  Schwinger  technique.)  Finally,  we  shall  give 
a  physical  interpretation  of  the  adjoint  equation  which  in  conjunction 
with  a  certain  symmetry  property  of  the  variational  forms,  will  lead  to 
a  reciprocity  theorem  whose  significance  for  the  ionosphere  will  be  ex- 
plained. 

To  sum  up,  therefore,  the  present  paper  contains  a  complete 
formulation  of  the  Schwinger  variational  procedure  as  adapted  to  the 
scattering  problem  in  question.  Actual  conrputatlons  of  the  matrix  . 
elements  for  specific  scattering  media  representing  the  earth's  iono- 
sphere are  to  appear  in  the  future.  .  ■   '~ 
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2*  Statement  of  Problem 


rig.  1 

Eef erring  to  Fig.  1,  0  ia  the  origin  of  a  right-handed  cartesian 
coordinate  system  (x,y,z)  with  the  positive  y  direction  into  the  paper.  The 
regions  z  <  0  and  z  >  a  are  assumed  to  "be  free  space  with  dielectric  con- 
stant e  and  magnetic  permeahility  m>  .  The  region  0  $^  z  ^  a  is  filled    •; 
with  some  material  medium  having  the  same  magnetic  permeahlllty  ^  hut 
characterized  "by  a  complex  dielectric  constant  which  most  in  general  he 
represented  hy  a  matrix,  c(z).  The  elements  of  this  matrix  are  known,  con- 
tinuous, complex  functions  of  z  which  depend  also  on  the  frequency,  al- 
though for  the  sake  of  brevity,  this  latter  dependence  will  not  he  explic-  : 
itly  indicated.   (In  the  special  case  where  e(z)  is  independent  of  fre- 
auency,  its  elements  are  all  real  functions.)  We  shall  further  require 
that  c(0)  s  c(a)  »  c  Z,  where  I  is  the  unit  matrix;  hence  both  the 

0 

dielectric  constant  and  magnetic  permeability  are  coatinaous  across  the 
boundary  planes  z  ■  0  and  z  =  a.  '  ■ 


*  The  MKS  syrtem  is  used  throughout. 
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We  now  suppose  that  aji  elliptlcally  polarized  plane  electromag- 
netic wave  with  harmonic  time  dependence  of  angular  frequency  w  impingeB 
on  the  medium  in  question  from  the  region  z  <  0.   In  Fig,  1,  the  paper  is 
the  plane  of  incidence,  L-L*  is  its  intersection  with  a  typical  incident 
plane  wavefront,  and  the  angle  of  incidence  is  Q,   measured  clockwise  from 
the  positive  z  direction.  We  seek  the  intensities  and  polarizations  of 
the  waves  reflected  hack  into  the  region  z  <  0  and  transmitted  into  the 
region  z  >  a. 

To  state  the  problem  more  completely:  1®*  ^H  field  vectors 
have  time  dependence  e~  *"  and  let  E  and  H  represent  the  spatially  de- 
pendent parts  of  the  total  electric  and  magnetic  fields  respectively  at  any 
point  (x,y,z), where  E  and  H  are  to  ohey  the  following  conditions: 

a)  E  and  H  shall  he  continuous  with  continuous  first  derivatives 
and  shall  satisfy  Maxwell's  equations  throughout  the  whole  space, i.e., 

1)  7  XE  -  liuuH  =  0 

2)  7x1+  iw€E=  0 

where  ^  »  ^  everywhere  and  the  matrix  e  varies  as  previously  described; 
^  and  c  are  thus  everywhere  continuous. 

b)  Tor  z  <  0,  E  and  H  shall  have  the  forms 

=•   ik(x  sinO+zcoseV  =■   ik(x8in©-zcosO)  ^  ,  =   lk(xsinO  +  zcosO).  =  „lk&C8in^zooBO) 
j!<  e  +  J!i_  e  and.  a    e  *o^ 

0  E  o  TT 

respectively,  where  E  .  E_,  H  ,  and  EL  are  constant  complex  vector  aoplitudes 

and  k  «  w\/p>  c  •  E  and  H  are  the  amplit-udes  of  the  incident  electric  and. 
^'^o  0     0      0  _     __ 

magnetic  fields  and  are  assximed  known;  E_  and  K.   are  the  corresponding  ampli- 
tudes of  the  reflected  fields  and  are  unknown. 

c)  Eor  z  >a,  E  and  H  shall  have  the  forms  I^   gilc(ialn»fzcos«)  ^^ 

IL,  e   ^       ^  w  /  respectively,  where  E_  and  H_  are  the  vector  amplitudes 

of  the  transmitted  electric  and  magnetic  fields.  Like  E_  and  H_,  E_  and  H_ 

are  unknown. 

Then  the  specific  problem  is  to  obtain  E^,  H^,  E^,  and  H^,  being 
given  the  vectors  E  and  H  and  the  fact  that  E  and  H  satisfy  conditions  (a), 

0      0 

(b),  and  (c).  It  should  be  noted,  however,  that  condition  (b)  is  easily 

*  In  eqxiation  (25t  c£  means  matrix  oniltiplication  of  c  and  £,  where  E  is  re- 
garded as  a  column  matrix. 
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shown  to  Imply  that  the  incident  and  reflected  electromagnetic  fields  in  the 
region  z  <  0  must  separately  satisfy  Maxvfell's  equations  and  that  therefore 
E  and  H  cannot  "be  assigned  independently  of  each  other  but  aust  "be   so 
chosen  that  the  incident  field  obeys  equations  (l)  and  (2)  in  the  region 
z  <  0. 

That  the  above  problem  is  well  set,  i.e.,  that  the  conditions  (*), 
(b),  and  (c)  are  Just  sufficient  to  determine  E  and  H  can  be  seen  in  a  some- 
what rough  way  as  follows:  from  Fig.  1  and  the  description  of  the  problem, 
it  is  clear  that  E  and  H  are  Independent  of  y.  Furthermore,  since  e  as  a 
function  of  position  depends  only  on  z,  the  assumption  that  E  and  H  have  an 
X  dependence  e       within  the  scattering  medium  as  well  as  outside  it 
transforms  eq-oations  (1)  and  (2)  into  a  system  of  ordinary  differential 
equations  in  z  within  the  scattering  ^ledivim.  By  simple  eliminations,  this 
system  can  be  reduced  to  a  set  of  four  linear  first  order  equations  with  con- 
tinuous coefficients  in  which  the  dependent  variables  are  the  x  and  y  com- 
ponents of  E  and  H.   It  is  well  Icnown  that  such  a  system  always  possesses  a 
continuously  differentiable  solution  in  which  each  of  the  four  unknown  func- 
tions depends  linearly  On  the  four  constants  of  integration  of  the  system; 
these  foar  constants  are  of  cotirse  to  be  determined. 

Now  ve  have  seen  that  outside  the  scattering  medium  the  reflected 

and  transmitted  fields  must  individually  satisfy  Maxwell's  equations,  so 

that  of  the  four  vectors  BL,  L,E_,  and  E_,  only  the  last  two,  for  example 

need  be  solved  for.  Moreover,  the  reflected  and  transmitted  fields  have  been 

(13)      _     _ 
assumed  to  be  plane  waves,  for  which  it  is  known     that  E_  and  E_  must  lie 

in  the  planes  of  their  respective  wavefronts.  Therefore  E_  and  E™  can  possess 
only  two  independent  components  each.  Thus  the  field  outside  the  scattering 
medium,  like  the  field  within  it,  is  a  linear  function  of  four  independent 
unknown  constants,  making  a  total  of  eight  ^Inknown  constants  to  be  determined 
in  the  problem.  Since  the  field  vectors  are  to  be  everywhere  continuous,  how- 
ever, imposition  of  the  continuity  requirement  on  the  x  and  y  components  of  E 
and  H  at  the  boundary  planes  z  =  0  and  z  =  a  leads  to  a  non-homogeneous  sys- 
tem of  eight  linear  algebraic  equations  for  the  eight  vinknowns.  The  deter- 
minant of  the  system  can  be  shown  not  to  vanish,  and  therefore  there  is  one 

and  only  one  solution, _______«_______ 

13.  Stratton,  J.  A.  Electromagnetic  Theory  p, 272-273. 

*  The  system  is  non-homogeneous  because  E  and  S  are  known.  " 

0      0 
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3.  Derivation  of  the  Integral  Equation 

In  the  present  section  we  shall  derlre  the  fundamental  rector  inte- 
gral equation  for  the  electric  field  In  the  scattering  region  0  5  z  *  a,  this 
eqiuition  constituting  the  actual  starting  point  of  the  theor7.  Specifically, 
we  shall  show  that  if  E  and  H  satisfy  conditions  (a),  (h),  and  (c)  of  Section 
2,  then  the  z  dependent  part  of  E  must  "be  a  solution  of  the  integral  equation. 

The  first  step  in  our  procedure  is  to  construct  a  Hertz  potential 
vector  to  which  ST  and  H  are  related  "by  a  pair  of  differentiation  formulas, 
and  to  this  end  we  employ  the  follw^ing  principle    in  electromagnetic 
theory!  let  the  vectors  f  and  g  "be  any  solution  to  Ma3well»e  equations  in 
a  region  where  ji  «  ^4  and  e  is  a  continuoua  "but  otherwise  arbitrary  matrix 
function  of  position.  I.e.,  , 

3)  7  X  f  -  iJUM-Qg  =0       ' 

k)  7  X  i  +  ituef  «  0 

Or'  X 

Then  there  always  exists  at  least  one  vector  function  9  satisfying  the  equa- 

'  "*  . 

tion 

2r     1 


5)       ^^   +k2f  =-  l-.[-e-c^I^i 


c 


and  having  the  farther  property  that 

6)  f  -  ^(^  •^)  +  k^^ 
and  _ 

7)  .        i  =  -i«.c^^^^ 

In  equation  (5),  7^  represents  th6  ordinary  laplacian  operating  - 
o6  each  of  the  rectangular  compSnents  of  f  ,  k  is  ui\/^Q*Q      as  "before,  and 
I  is  the  unit  matrix.  0  is  called  the  Hertz  potential  of  the  electromagnetic 

field  defined  hy  f  and  g.  »  >' 

The  foregoing  principle  evidently  assures  us  that  corresponding  to 
the  pair  of  vectors  B(x,z)  and  H(x,z)  obeying  conditions  (a),  (b).  and  (c)  o£. 
Section  2,  there  always  exists  at  least  one  function  if(x,«)  say,  satisfying 
the  equation 

8)  ^  TU.z)  +  k^  ?(x.z)  -  -  ~  C«(0  -  V^  ^^**'^     ' 


•  This  principle  is  proved  in  Appendix  I. 
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for  all  X  and  t,  emd  to  which  E  and  H  are  further  related  "bj   the  formulas 

9)  I(x,?)  -  7  [>  •  v(3c.>ra  +  k^  ijrCx.O 

10)  H(x,»)  -   -lu.<o  7Xv(i.e) 
for  all  X  and  z. 

Ve  shall  now  show  that  the  existence  of  a  function  y-  which  satis- 
fles  e(i'uation8  (8),  (9)»  and  (10)  necessarily  leads  to  the  detemioation 
of  its  form,  aad  that  from  the  form  of  ^  we  can  deduce  that  the  t   dependent 
part  of  X(x,p)  satisfies  an  integral  equation. 

To  begin  with,  consider  the  function  ^  outside  the  soatterlog 
region.  There,  the  right  side  of  eq\:iation  (8)  mast  Tanlsh  since  «(z)  =  'q-*^ 
and  hence  conditions  ("b)  and  (c)  on  E(x,z)  and  E(x,2)  together  with  equa- 
tione  (9)  and  10)  inrply  that 

,,v  —,.       -       ik(xBinO  +  zco80)  .  -  ^ik(x8in« -jtcosfl)        ^- 

11)  yf{x,%}   ■  v^  e  ^  '''e  '      z  <  0 

,_^   -/   V   -   ik(x8in«  +  zco80)  ^„ 

12)  v(x.z)  =  Vj  e  '  z  >a 

where  \f  ,  \j?-,  and  vjr-  are  consteuit  complex  vector  amplitudes. 

0    A        T 

With  the  aid  of  equations  (11)  and  (12),  which  express  the  form  of 
^  outside  the  scattering  region,  we  can,  in  the  followisg  manner,  obtain  an 

expression  for  \^  which  is  valid  anywhere » 

-     ~"     ikxsind 
Semembering  that  the  x  dependence  of  £  and  H  is  e       ,  it 

follows  that  \^  must  have  the  same  x  dependence.  If  in  eqixation  (8)  we  denote 

the  vector  HeCz)  -  c  ll  E(x,z)  by  P(x,z),  the  introduction  of  the  x  depend- 

^  0   "^ 

ikxsinO   ^     j.^      i.j   *. 
ence  e       reduces  the  equation  to  ,  ,  .    .,  .,. 

"2 

13)  ^-^+  k^  cob2«v(z)  =  -  i-P(z) 

dz  0 

where  v(z)  and  P(z)  are  the  z  dependent  parts  of  ^(x,z)  and  P(x,z), 

Equation  (13)  of  course  represents  three  separate  scalar  equations, 
each  of  which  involves  only  one  of  the  three  components  of  ^(z)  and  P(z). 

Therefore,  we  may  Introduce  a  scalar  Greenes  function   0  of  the  form 
...  ikiz-z'icos  Q 

;^ '"'•->-  °  ziL. 

1^.  Courant-Hilbert,  Method,  der  Hath.  PhyB.I,p.302.       -f::   ;,■  :■=    : 

■■■■■■■■     /  ■  ■-  "  '  ~  "'      ■■-'•-  -— —  •■  _ 
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and  apply  Green's  formala   to  v(z)  and  G(z,z')  over  the  interval  -oo  ^  z  ^  oo. 

There  results,  after  changing  the  Independent  ▼aria'ble  from  z  to  z' 
in  equation  (13). 


CO  / 

15)'^    r   G(z.z«)P(z')dz«-l^^G(z.z.)-^2^f^i(z«) 


'-00 


Z'=  2-0 


-00 


z'=oo 


z'=E+0 

At  z'^-oo,  the  form  of  vp(z')  is  evidently  given  hy  the  z  dependent 

ik(z-8')  C08« 


Insertion 


part  of  eapression  (11)  while  the  form  of  G  is  2ik  cosO 

of  these  expressions  into  the  first  term  on  the  right  of  eq:aation  (15)  gives 
as  the  lower  limit  of  that  term  ~ 

ik(s-zi)oo8  0 

— XK.Z  "  <JOB  Vl  «' 


..K        r,  «/-     ikz'  cosO      -  ^-ikz 

16)  -     ik  cosOl\|r  e  -  Vx,e 


2  ik  cosO 


^  1     lk(z-z»)cos«/-  ^ikz» 
*2^  VV 


008  ©  ^  -     -ikz 


l008  0\ 


which  reduces  simply  to 


17) 


-       ikz  cosO 
—  w     e 
*o 


Similarly,  at  z'=  +oo,  \y(z')  "behaves  like  the  z  dependent  part  of 

ik(z '-2)008  0 

expression  (12)  while  G  takes  the  form  2  ikoosO —  *   Insertion  of  these 

expressions  into  the  second  term  on  the  right  of  equation  (15)  gives  as  the 
upper  limit  of  that  term 


fQ\  4,     ft-  ikz •  DOS 0  e^^^'  -zjeoa  «   ^  ^^^^^   _z)coe  0_  ^ikz  cos  0  _  ' 

18)  Ikcose^ie  2  ik  cos  9  '  Z  ^  ^T  ^         =  °- 

Dae  to  the  continuity  of  G(z,z')  at  z'=z,  the  terms  on  the  right  of 
eqjiation  (15)  involving  t^  must  cancel  at  z'=z,  leaving  only 

z«=a-0 

19)  .S^J^^(„)  I  '„.  [.ie"^(^-='>°°«%(z.)-|e"^(^'-«)'=°«V')]  ""^^^^ 


z'»24-0 


Z'=Z 


15.  Ihid. ,  p.  237. 


-c 


(:'■ 
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Thus  equation  (15)  Day  now  "be  vritten 


It  follows  from  its  definition,  howerer,  that  the  vector  P  vanlshe* _ 
outside  the  region  0  ^  s  ^  a.  Therefore,  ths  integration  in  equation  (20) 
need  extend  only  from  0  to  a  so  that  the  equation  hecomes 

a   I    I 

olN      -/    \         -   ^ifczCDSO  1 /■       Ik    Z-Z«    COS©     =f     ,v     .     , 

a)     v(z)=,p^e  -     2lckoosfiy    «  '  P(^')dz' 

0       0 

where  we  have  replaced  G  hy  its  explicit  form,  given  by  eq\as.tion  (1^)» 

Moltipllcation  of  hoth  sides  of  equation  (21)  by  e       results  in 

o*\      -/   \   -  ik(xalnO+  zoo8«)    e^^^°^^^   /^ik|z-z«  IcosO  ^/  ,v  .  , 

22)  ,,,(x.z)-Vo«  -  2ic  kcDsOy  ^  P(z')dz» 

0       0 

One  further  necessary  condition  on  the  form  of  y(x,z)  can  he  deduced 
from  condition  (b)  of  Section  2  and  equations  (9)  and  (11).   It  is  that  ijr  , 
which  we  have  not  yet  specified,  must  be  such  that 

23)  .      [^(V  .  )+  k^l  y^  ,ik(iBln«+zoo8«)  „  |^  ^ik(xsinO+  zcob«) 

A  breakdown  of  equation  (23)  into  its  scalar  components  together 
with  the  fact  that  E  must  lie  in  the  plane  of  the  incident  wavefront  shows 
at  once  that  w  can  always  be  chosen  so  as  to  satisfy  equation  (23) •  and 

0 

indeed  in  an  infinite  number  of  ways. 

Our  results  at  this  point  inay  be  3\iinmed  up  as  follows:   if  E(x,z) 
and  H(x,z)  obey  conditions  (a),  (b),  and  (c)  of  Section  2,  then  the  vector 
^(x,z)  telated  to  B(x,z)  and  H(x,z)  through  equations  (8),  (9).  and  (10) 
must  necessarily  be  of  the  form  given  by  equation  (22)  for  all  x  and  z^ln 
which  ^  is  any  vector  satisfying  equation  (23). 

Ve  may  now  proceed  directly  to  the  deduction  of  the  Integral  equa- 
tion as  a  necessary  consequence  of  the  restilt  Just  stated.  Applying  the 
operator  lv(7  >)  ^  Ic  j  to  both  sides  of  equation  (22),  we  obtain  in  view 
of  equations  (9)  and  (23) 


»)  i(x,.) .  v<^''-'«*  '"■*'  -  [»(»•  >  ^'izirto'iA'^'""'"'^"'  *" 
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where  the  effect  of  the  operator  on  the  second  term  on  the  right  remains 
to  be  worked  out.  This  calculation,  although  elementary,  is  rather  com- 
plicated end  we  have  therefore  relegated  it  to  Appendix  ZI.  The  final 
result  is 

0       0 

a 


^j^^ikxaln© 


/^Iklz-z'lcos©  .  /   ,v-,  ,x.  .  .  1  .ikxslnO^  /  v^. 


o 
in  which  A(s,zi)  is  a  matrix  defined  "bj 


2e    oosD        . 
o  -'.  .  o 


2 

cos   0  0  -slnOoosO 


26)  A(z,B«)  =1  0  1  0|««<z 

,-8in  9  cos  9-  0  sin     Q 

COB  Q  0  sinOcosd  \ 

27)  A(z,z«)  =1         0      1     .  0  I  ze  >  z 

BiirO  COS  0      0  sin  G  / 

end  in  which  AP  in  the  integral  means  matrix  molt ipiicat ion  of  A  end  P  regerd- 
ing  the  vector  P  es  e  column  matrix,  P^(z)  is  the  z  component  of  P(z),  and  k 
is  a  omit  vector  in  the  positive  z  direction. 

Suhstituting  for  the  second  term  on  the  right  of  equation(2^)  the 
evaluation  as  given  by  equation  (25).  we  get  after  dropping  the  x  dependence 


28)  i(z)  ^  i-P3(z)  k  =  E^  ,ikz«B«^  __^y%lkh-«'|cos«^(^^^,jj(^,j^^, 

0 

It  is  convenient  to  introduce  two  new  matrices,  B(z)  and  C(z).  B(z) 
is  defined  through  the  relationship 

29)  CqBCz)  =  €(z)  -  e^I 

and  therefore  B(z)  vanishes  outside  the  region  0  ^  x  ^  a.  From  the  definition 
of  P(x,z)  we  see  that 

30)  P(z)  «  e^B(z)  5(z) 

B(z)  is  the  so-called  electric  susceptibility  matrix  of  the  scatter- 
ing medium,  while  P  la  known  as  the  electric  polarization  vector  of  the  medium. 

The  matrix  C(z)  Is  defined  as  having  its  first  two  rows  zero  and  its 
third  row  eqtial  to  the  third  row  of  B(z).  Evidently  then        .  ,  ,.,  .* 

31)  '   P3(z)  k  =  e^CCz)  E(z) 
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Substituting  eq.tiatlona   (30)  and   (3])      Into  equation   (28),  ve 
finally  obtain 

a 

32)  Cl*  C(z)J  i(0  =i/^'=«'««*  j^/e^^l^-'^'l^''««A(z,z«)3(..)E(z')d.. 

0 

Equation  (32)  is  the  desired  rector  integral  equation  for  the  z  depend- 
ent part  of  the  electric  field  in  the  scattering  region  0  ^  2  ^  a.  Together 
with  the  adjoint  integral  equation  referred  to  in  the  Introduction  and  to  be 
derived  presently,  it  plays  a  fundamental  role  in  the  construction  of  the  vari- 
ational formulas  for  the  elements  of  the  reflection  and  transmission  matrices. 

There  axe  several  points  regarding  equation  (32)  which  we  now  wish 
to  mention  briefly.  First,  although  the  integration  extends  only  over  the 
scattering  region,  the  method  of  deriving  the  equation  implies  that  it  cor- 
rectly represents  the  electric  field  anywhere;   in  fact,  the  equation  has  the 
clear  physical  interpretation  that  it  expresses  the  total  electric  field  at 
any  point  as  the  sum  of  the  incident  field  and  a  field  produced  by  the  inte- 
grated effect  of  induced  soiarces  within  the  scattering  medium. 

Second,  equation  (32)  is  a  system  of  three  scalar  integral  equationa 

for  the  x,y,  and  z  components  of  l!(z)  and  as  such  it  possesses  a  well  estab- 

T6  17 
liahed  existence  theory  '   .  This  theory  assures  ua  that  the  equation  has  a 

unique,  continuously  differentiable  solution,  because  the  kernel  is  continu- 
ous except  for  th6  finite  jump  in  A(z,a')  at  z  =  z'.  Now  we  have  Just  fin- 
ished showing  that  every  electric  field  obeying  tho  conditions  of  our  original 
problem  provides  a  solution  to  equation  (32),  and  the  uniqueness  theorem  for 
the  Integral  equation  now  implies  that  the  converse  statement  is  also  true, 
viz.,  if  we  start  with  equation  (32)  in  which  the  first  term  on  the  right  I0 
the  z  dependent  part  of  a  plane  wave  solution  to  the  free  space  MEiiwell'o 
equations,  then  the  solution  to  the  eq'uation  fximishes  an  electric  field        , 
satisfying  all  the  conditions  of  our  original  problem.  Thus  the  integral 
equation  and  the  original  problem  are  fiilly  equivalent, 

rinally,  it  should  be  noted  that  the  matrix  A(b,zI)  is  singular; 
the  significance  of  this  fact  will  become  more  apparent  in  the  next  section, 
dealing  with  the  reflection  and  transmission  matrices,  where  we  shall  deduce 
several  iaiportant  and  useful  properties  of  A(z,z')  associated  with  its  singu- 
lar character, 

16)  Helllnger  and  Toepllti,  Bncyklop.Math.Wiss.  Bd.  2,  Tell  3,  Heft  9. 

17)  Platrler,  Ch.  Journal  de  Math.  Series  6,  vol. 9.  p. 257  ff.  1913. 
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k.   The  Hef lection  and  TransmisBlon  Matrices 

In  the  Introduction,  we  allnded  to  the  fact  that  the  prohlem  of  o"b- 
taining  the  complex  vector  amplitudes  E_  and  E_  (and  thereby  automatically 
Hp  and  HL  )  can  "be  reduced  to  that  of  calculating  the  elements  of  a  certain 
•reflection  matrix*  and  a  certain  •transmission  matrix"  which  operate  on 
B  to  yield  E^  and  K,  respectively.   In  this  section  we  shall  carry  out  the 
reduction  and  relate  these  matrices  to  the  integral  equation  (32). 

We  "begin  "by  considering:  any  two  unit  orthogonal  vectors,  u-  and  u-. 
In  the  plane  of  the  incident  wavefront.   Observing  once  again  that  for  electro- 
magnetic plane  wav^s  in  free  space,  the  electric  and  nBgnetic  vectors  lie  in 
the  plane  of  the  wavefront,  we  may  evidently  write  >- 

33)        \'Ci^^h*^i'^^ 

where  t^  '  and  (^   are  complex  scalare. 

Similarly,  if  u,  and  u-  are  any  two  unit  orthogonal  vectors  in  the 
plane  of  the  reflected  wavefront,  we  have 

3.)  i,=t(^>ui*e<^'u'  - 

where  C^  'and  C  ^  'are  complex  scalars, 

Finally,  since  the  planes  of  the  transmitted  and  incident  wavefronts 
are  parallel,  it  is  convenient  to  choose  the  vectors  u.  and  u.  of  equation  (33) 
as  unit  orthogonal  vectors  in  the  transmitted  wavefront.  There  follows 

35)  5,  =  £(^)5,*ef :, 

where  of  course  ^  J,     and  C^  'are  complex  scalars. 

The  original  pro"blem  posed  "by  conditions  (a),  ("b),  and  (c)  of  Section 
2  is  linear;  hence  the  effect  of  multiplying  i  "by  a  complex  n'um"ber  is  simply 
to  multiply  R-  and  E_  "by  the  same  num"ber,  while  the  effect  of  superimposing 
two  incident  waves  in  the  original  prohiem  is  to  superimpose  the  corresponding 
reflected  and  trensraltted  waves.  These  consequences  of  linearity  then  lead  us 
immediately  to  the  following  relationships: 


36)       : 

c  (1)  .p      t^^KiL     C    (2)               ' 

37)      ■'' 

■^-    --   C      (2)    _   -B          C   (1)    a.   fl          C      (2)         -i;    a  !-'<." 

*=E       -  ^21  *=-o       *  ^22  Co 

38)    ■■; 

■  r:    C  T       "  ^11  ^0       *  ^12    Co             :;-?  ,j    r 

39) 

g   ^)=T       ^^^)  +  T        C    (2) 
^T           ^21  '-■0      *  ^22  *=-o 

•  T 
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where  the  R. .  and  T.  .  are  complex  scalars.   In  matrix  form,  equations  (36) 


and  (37)  may  "be  written  as 

40)  Eg  =  R  i^ 

and  equations  (38)  and  (39)  as 
hi)  ij  =  T  i^ 


The  matrix  R  = 


|r   II  is  the  so-called  reflection  matrix  of  the 


scattering  medium,  while  T  b||;  ||  is  the  transmission  matrix  of  the  medium. 
In  addition  to  depending  on  the  specific  nature  of  the  scattering  msdium  as 
characterized  "by  its  susceptibility  matrix  B(z),  the  R. .  and  T .  .  are  func- 
tions of  the  angle  of  incidence,  Q,  s^d  the  frequency,  m »    Obviously,  a 
knowledge  of  the  elements  R. .  and  T.  .  constitutes  the  solution  to  our  prob- 
lem;  thus,  the  reflection  and  transmission  matrices  contain  all  the  desired 
information  relevant  to  the  effect  of  the  scattering  meditun  on  the  Incident 
wave. 

How  our  eventual  object  with  regard  to  the  elements  R.  .  and  T .  . 

is  the  derivation  of  a  variational  principle  for  them,  but  in  order  to 
accomplish  this,  we  must  first  relate  them  to  certain  particular  solutions 
of  the  integral  equation  which  correspond  to  special  forms  of  the  incident 
wave.  We  begin  by  examining  the  integral  eqioation  (32)  for  z  <  0  and  z  >  a,  i.e., 
outside  the  scattering  region.  Taking  first  the  case  of  z  <  0,  we  see  that 
C(z)  =  0  and  |z-z'|  =  (z»-z),  while  A(z,z')  takes  the  form  of  the  constant 
matrix  given  in  equation  (27).  which  we  denote  by  A  .  Equation  (32)  then 
becomes 

a 

ii«\  Zf   \       ^       ikzoos©  .     -ikzcos^  ik         Z'  ikz'oosO  .  _/    ,%->■    ,\,    ,         ^  « 

'*2}  E(z)  =  E^  e  ■>■  e  ^  ^^^^    /  e  A_^B(z«)E(8' )dz«,   z  <  0 

0 

•  From  the  form  of  the  rig^t  side  of  equation  (42)  together  with  the  con- 
sideration of  condition  (b)  of  Section  2,  it  is  clear  that  we  may  write 

0 

^or  the  case  z  >a,  we  again  have  C(z)  =  0,  while   |z-2'|   »   (z-z') 
and  A(z,z»)  takes  the  form  of  the  constant  matrix  given  in  equation  (26), 
which  we  denote  by  A.  «     Equation  (32)  then  becomes 

iii.\  5/  \      ^      IkzoosO  .      Ik  zoos  ^      ik       r  ^-iks'cosO  .  -^   ,*-#   ,».    ,        ,. 

*W)    I         E(z)  »  E^  e  +e  o^,a    /  ^  A  B(z«)E(t')dz«,   s  >  a 

r  I   .  O  2008  9    /  ~ 

•      •;■     ■  ...    ^.  'O      • 
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and  from  the  fora  of  the  right  side  of  eq.uation  (W)  together  with  condition 
(c)  of  Section  2,  we  get 

0 

Tormulas  (^3)  ^^^d  (^5)  expresa  the  vector  amplitudes  of  the  reflected 
and  transmitted  waves  in  terms  of  the  electric  field  in  the  scattering  region 
0  <  z  <  a. 

How  let  E/.»(z),  (j  «»  1»2),  represent  those  solutions  of  the  inte- 
gral equation  (32)  which  correspond  respectively  to  the  incident  waves  of 
vector  amplitiade  1  /.v  =  n.,    (J  «=  1,2),  and  let  S_y  .v  and  ^/j\  represent 

the  corresponding  reflected  and  transmitted  wave  eunplitudes.  Then  accord- 
ing to  equation  (43) 

0 

and  according  to  equation  (45) 


0 

From  equation  (33)  we  have 

where  from  the  definition  of  E  *.v  it  is  evident  that 

ixQ^   '  ?1^^   „  .   _  1.  ?c=  J  (k-1.2) 

^^  ^^J)  °  \j  -  0.  k?<  J  (J=l,2) 

Substitution  of  this  last  relation  into  equations  (36)  and  (37)  gives 


(J)    k=l  ^^     ^^         ^^  (J«l,2) 


50)         £.^;   =  2:  Rik  8kj  ■=  H^j 


omee 


and  since  equation  (3^),  when  E_  =  E_/.\.  "bee 

it  follows  from  equation  (50)  and  the  orthogonality  of  the  vectors  xu   and  yx-   that 
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Slinllarl7,    if  we  su^bstitnte  equation  (^)  into  equations   (38)  and 
(39),  we  get 


(1)  2  (1»1.2) 

»  >      T       6       »  T 

(J)  fel  ^-^  ^J  ^J  (J- 1.2) 

and  since  equation  (35),  when  E_  =  Il_/^,  becomeB 


53)  £^^  »  T  T   6   »  T 


^^)       ^o)=^^::.^i^^^i- 


(J"  1.2) 


(J)  "     "(J)  ^ 
it  follows  from  equation  (53)  <Mid  the  orthoeonality  of  u-  and  u, 

'•'^         ^  (J- 1.2) 

Bqijations  (52)  and  (55)   are  significant  'because  they  express  the 
essential  meaning  of  the  reflection  and  transmission  matrix  elements.  Thus 
we  see  that  H. .  represents  the  component  in  the  direction  u.  of  the  vector 
amplitude  of  that  reflected  wave  which  was  excited  by  an  incident  wave  of 
vector  amplitude  u.{  naturally,  T^.  has,  mutatis  mutandis,  the  same  meaning. 

g 

If  we  now  take  the  scalar  product  of  equation  {h6)   with  u.  ,  we 
obtain,  in  view  of  equation  (52) 

o 
Likewise,  the  scalar  product  of  equation  (^7)  with  u.  produces.  In 
conjunction  with  equation  (55) 

a  (1=1.2) 

^^>      ''ij  '  «ij  *  ^^«"'''''*'''  ^  •  ^-  ^(-'^  ^J)(^')  ^*»  (J=1.2) 

It  is  possible  to  effect  a  useful  simplification  of  formulas  {56)   and 
(57)  by  means  of  the  following  property  of  the  matrices  A  and  A  :  Every  vec- 
tor  in  the  plane  of  the  reflected  wavefront  is  an  eigenvector  of  A  with  eigen- 
value  1,  while  every  vector  in  the  plane  of  the  incident  or  transmitted  wave- 
front  is  an  eigenvector  of  A_  with  eigenvalue  1» 

To  prove  this  statement ^  let  P  be  any  vector  in  the  plane  of  the 
feflected  wavefront.  Then  in  terms  of  its  rectangular  components,  7   most  / 
clearly  be  of  the  form 

58)  ^        _   f  -  r,  coB«  i  +  ?  > y  sin^  k    '    '.''■  ' ' 


!"■  ■  1  '. 
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where  F,  and  F-  are  arbitrary  complex  sealare.  The  product  A^F  glres 


59) 


V  = 


cos  0  1     cos«+  sln^  coB«r.    sin^t     /  I^cobQ 

7r. 


BinQ  coaQT.   coaQ*  sin  OT.slnO 


=    F, 


■  F 


F^  8ln«i 


which  proves  the  statement  for  A  . 

Similarly,  let  0  "be  any  vector  in  the  incident  wavefront.  Then  C 
may  he  written 

60)  G  «  Oj^  co8«  i  +  Ggi  "  ^1  8iJ**k 

where  G,  and  G-  are  arbitrary  complex  scalars,  while  the  product  A_G  hpcomes 


COB  0     0      -Bin 4  cos  4 


jO^coe  q\ 


61) 


A  G 


0     1 
-sin4  cos  Q     0 


Bln^  0 


l-Gj^eln  0 


'/ 


coB^O  6^co8«-  sln«co8«  (-Gj^sln«)\     (q^  cobO 


=    G, 


k-sinfi  co8«G^cos«+  sin  0(-G^8ln«)    /     l-G^8in« 


Thus,   the  statement  is  proved.     We  apply  this  result  hy  noting  first 
that  it  implies 

62) 


A^u^  -^ 


(i=1.2) 


63) 

once  that 
6h) 

65) 


K  h    =  "^i 


(1=1,2) 


Now  A  and  A  are  "both  symmetric  matrices,  from  which  it  follows  at 
u|  •  A^B(z»)  B^^j(z«)  «=  A^.  u^  .  B(z»)  E(j)(2»)  -;;  > 
u^  .  A_B(b»)  E^  j(z«)  -  A_  u^-  B(z»)  i^jj(z«)   ^°^^ 
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whence  from  equations  (62)  and  (63)  we  obtain 


«>  •   hi  ■  j^.A^'"''  <  ■  =('•>  ^.)(-5 '" 


(J=l,2) 


0 

a 


o 
Equations  (66)  and  (6?)  constitute  the  desired  form  of  the  relation- 
ehip,  already  referred  to,  between  the  particular  solutions  E/.v(z)  of  the 
integral  equation  (32)  and  the  elements  of  the  reflection  and  transmission 
matrices.  Although  these  relations  may  themselves  be  used  for  approximation 
of  the  E   and  T  .  by  replacement  of  E/.v(x)  in  the  integrals  with  suitable 

approximate  trial  functions,  their  chief  importance  for  ub  will  be  in  the 
fundamental  role  they  play  in  the  derivation  of  the  variational  formulas. 

Before  proceeding  with  the  derivation,  to  be  carried  out  in  the 
next  section,  we  wish  to  interpret  the  singular  character  of  the  matrix 
A(z,2').  I'or  this  purpose  we  need  a  farther  property  of  the  matrices  A^ 
and  A  ,  viz.,  every  vector  perpendicular  to  the  reflected  wavefront  is  an 
eigenvector  of  A  with  eigenvalue  0, while  the  corresponding  statement 
holds  for  A  and  the  incident  wavefront. 

To  verify  this  assertion,  let  V  be  any  vector  perpendicular  to  the 
reflected  wavefront.  We  may  write  7  as  ~~ 

68)  ?  =  V  sln^  1  -  7  cos«k 

where  7  is  any  complex  scalar.     Then 

9  \  /  \         /        2  J.. 

COS  «     0     sin«  cosOXf  7  8in«\       (cos    «7  8in^-8infico8«7co8« 


^sin^  COS  Q     0  sin   0/\-7  cos  fi/     \sin«co8«7  sind-sin  d  7coB^j   "^0/ 

and  in  en  exactly  analogous  way  we  can  show  that  -'p        ;,    U- 

70)  AW  =0  "'^-  ';;  '■    '  .  ' 

where  W  is  any  vector  orthogonal  to  the  incident  wavefront,  ,     '  ,'  ' 


aA  ir^.it-'to'-J  "Ire 


> 
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Now  recalling  that  A(i:,z»)  «=  A^  for  t  <  t*   and  A(8,z«)  =  A^  for  z  >  sS 

we  see  that  the  relationships  expresBcd  "by  eqriations  (59),  (61), (69)  and  (70) 
can  he  comhined  into  the  following  single  statement:   the  matrix  A(z,tO  pro- 
Jecte  erery  vector  orthogonally  onto  the  planes  of  the  reflected  or  incident 
wavefront  according  as  z  <  z'  or  z  >  z',  respectively.  Since  A(z,zO  thereby 

maps  a  three  dimensional  space  onto  a  two  dimensional  one,  it  must,  hy  a  well 

18 
known  theorem   in  linear  algehra,  he  singular. 

This  property  of  A(z,z')  enahles  us  to  verify  immediately  that  the 
vectors  IL  ani.  E_  as  given  by  formulas  (^3)  and  (^5)  actually  lie  in  the  planes 
of  their  respective  wavefronts  as  required  by  Maxwell's  equations;  the  geo- 
metrical significance  of  the  singular  character  of  the  matrix  A(z,z')  is  thus 
made  clear. 

To  conclude  this  section,  we  introduce  a  matrix  S  related  to  T  by 
the  equation 
71)  T  =  I  -f  S 

Evidently  then,  according  to  equation  (67) 

0 

It  will  prove  more  netural  to  obtain  the  variational  fomulas  for 

the  elements  S .  .   than  for  T.  .  .., 

IJ  Ij    ,..•'  ;v;    ^   ::,■'    ■''-■■■•    ■.■^'      \>  -i:  rip.  •...:'  ^--    .  ,.    ,  . 


«  I?         ..  t-  - 


'■m 


'^i'/".    ■■ 


< 


; 

.-  ■^n;^^ 

V 

'.iC^'- 

\      • 

. .:    6 

i      T.  i 

1  t    i 

- 

if       V     ™y  »•■:, 

.;■  ■'": 

■•  i  ■. 

•'X 


i..' ^ -.';"•  p  .-       ■.'•:■■■??•' .■;!,  .    \y    -.'i       •  :    -^    4.^  ft-   ^      >rvtf     .■.'•■'    ■;?   '-'  !ii'-/'i 


18*.  Birkhoff  and  Maclane,  A  Survey  of  Modem  Algebra,  p.   209i 
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5.  The  Variational  yonnalaa 

We  are  finally  in  a  position  to  derire  the  Tariational  fomralaa 
for  the  elements  of  the  matrices  E  and  S.  The  starting  point  of  the  deriv- 
ation is  the  construction  of  a  certain  variational  form  involving  inner 
products  of  elements  in  an  abstract  linear  space;  we  shall  prove  that 
this  form  is  stationary  with  respect  to  small  variations  in  the  elements 
on  which  it  depends  if  these  variations  take  place  in  the  neighborhood 
of  two  particular  elements  of  the  space.  Then,  by  constructing  an  ex- 
plicit realization  of  the  abstract  space,  we  shall  set  up  a  correspond- 
aace  between  the  elements  of  the  latter  and  the  vector  functions  of  our 
original  problem.   Insertion  of  these  vector  functions  into  the  variation- 
al form  will  lead  to  the  desired  stationary  formulas  for  R.  and  S. , 

Tor  convenience,  we  recall  briefly  the  definitions  of  linear 

19  20 
space,  inner  product,  and  linear  operation  '   .   In  all  of  our  considera- 
tions of  the  abstract  space,  smell  Greek  letters  denote  elements  of  the 
space,  small  Roman  letters  denote  ordinary  complex  numbers,  and  capital 
Roman  letters  denote  linear  operations, 

A  "linear  space  ■  is  a  set  of  elements  with  the  following  proper- 
ties: 

1)  Tor  every  pair  of  elements  C,  and  {_  °^  *^®  space,  there  exists  a 
commutative  end  associative  operation  indicated  by  +  such  that  t,  +  Cg  is 
also  an  element  of  the  space, 

2)  Corresponding  to  any  complex  number  c  and  element  t,  there  is  a 
distributive  and  associrtive  operation  denoted  by  ct,  such  that  ct  belongs 
to  the  space,  and  1«C  =  t.  • 

3)  There  exists  a  null  element  in  the  space,  denoted  by  0,  with  the 
properties  C  +  0  =  0,  c«0  =0,  O-t  =  0  ___ 

The  "inner  product"  of  two  elements  in  the  linear  space,  denoted 
"by  (CifCp),  is  a  complex  nunber  having  the  properties  that  (cCj^^.tg)**  "^  ^^p^a^* 

for  {.,  5^  0,  (C.  ,t-,)  =  0  for  {,  »  0.  The  asterisk  denotes  complex  conjugate. 

A  'linear  operator"  ,  K,  in  the  linear  space  is  a  correspondence 
which  maps  each  element  C  of  the  space  on  an  element  Kt  of  the  space  in  such 
fashion  that  for  all  elements  C  and  1   and  all  coinplex  numbers  c  and  d, 
K^cC  +  dl  )  =  cKC  •••  dZt  . 

19.  Birkhoff  and  Maclane,  Op.cit.,  Chapter  7. 

20.  Stone,  M.H.,  Linear  Trans,  in  Hilbert  Space,  Chapters  1  and  2, 
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We  define  the  linear  operator  which  is  %.djolnt"  to  K  and  which 
we  denote  hy  K',  through  the  property  that  for  all  elements  t  and  1.  , 
(K'C.  t)  =  (C.Kt  ). 

Finally,  if  L  and  K  are  two  linear  operators,  then  by  the  operator 
IK  applied  to  an  element  (  we  mean  L(EC}.  ^he  adjoint  of  IX.   ia  then  easily 
shown  to  "be  K'L', 

Now  let  C  and  C '  "be  any  two  variable  elements  of  the  linear  space, 

let  a  and  a'  he  two  fixed  elements  of  the  space,  let  L  and  K  be  linear  oper- 

21 
ators  and  L'  and  K*  their  respective  adjoints.   Consider  the  form 

73)    hrt£'l=  {T^n\o){0'l,U) 
n)        nLt.tJ-    (C'.LKLC) 

Trom  its  definition,  h  is  evidently  a  complex  number.  We  shall  show  that 
h  is  stationary  with  respect  to  arbitrary  small  variations  in  C  and  {' 
about  the  elements  oj  and  T\  ',  where 'Tl  satisfies  the  equation 

7^)  a  =  KIT) 

and  oi '  satisfies  the  equation 

75)  a'  =  E'l'^T)! 

More  precisely,  let 'h  andi}*  satisfy  equations  (7^+)  and  (75). 
Then,  if  for  t  and  C*  we  write 

76)  C  ='»)  +8f  C  »/»?'  +  B»f'  V, 

where  s  and  s'  are  complex  numbers,  and  {  and  f  •  arbitrary  elements  of  the 
linear  space,  we  shall  prove 

^^^       §7  ^  P*}*^.^  »^'  *  ^^ '1    "  ^  for  all  B 

8»=0 

78)       1^,  h  U+sJ  ,0J«  +  s»^ '"j   =0      for  all  8  8 

8=0 

Equations  (77)  and  (78)  express  a  property  of  h  which  is  stronger 
than  the  usual  stationary  one.  Normally  equation  (77)  would  hold  only  for 
8  »  0  and  equation  (78)  only  for  s'  »  0. 

To  prove  equation  (77),  we  first  note  that  from  the  defialtion 
of  the  adjoint  operator,  it  follows  that  (C'.IKLC)  »  (L«K»L»C«,t).  Intro- 

21,  lAx,  H.,  Phys.  Rev.  vol.  78,  p.  306,  May  1,  1950. 
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duclng  thlB  eqtiality  into  equation  (73)  and  setting  b«  «0,  we  get 

79)  h   [^+8|    ,/T)'J  =  (LiK'L'y.'>>+  sF  )  • 

Bat  according  to  equation  (75),  K'L'/);'  =  a'  bo  that  L'K'L'/i^'  =  L'a'. 

Therefore , 

80)  (L'K'L'O^  »,'r|+  8  E    )  =   (L'a'./^+Bn  "   (a.',  K'JJ+af  )) 

Stibstltution  of  the  last  term  on  the  right  of  equation  (80)  for 
the  denominator  in  equation   (79)  gives 

81)  h  [o;j+8^,^«]  »  (L'o^'.a) 

and  since  the  right  aide  of  equation  (8l)  ie  independent  of  a,  equation  (77) 
is  thereby  proved. 

The  verification  of  equation  (78)  is  accomplished  in  similar 
fashion.  We  set  s  =  0  in  equation  (73)  to  ohtain 

Using  equation  (7^),  the  denominator  in  equation  (82)  is  seen  to  obey  the 

equations 

83)      (/9«  +  B«  fe  »,  IKL-J))  «  ('9  •  +  8«|  »,  La)  »  (L'CP'  +s'i  «).  a) 

from  which  equation  (82)  "becomes 

Qk)  h[o;.'»?«  +  B«^«]  =  (a»,L'9) 

which  verifies  equation  (78). 

The  value  of  h[t,C']  at  the  stationary  point  C  ='>}.  C'  =^'.  oay 
he  written  in  several  forms  which  we  set  down  for  future  reference 

The  next  step  in  the  direction  of  applying  the  foregoing  stationary 
property  of  h  is  to  construct  an  ejpllcit  representation  of  the  abstract 
linear  space.  For  the  representors  we  shall  choose  the  set  of  all  three 
component  complex  rector  functions  of  the  real  variahle  z,  continuous  over 
the  interval  0  ^  z  5  a  .  Since  the  z  dependent  parts  of  the  electromag- 
netic field  vectors  in  our  original  problem  belong  to  this  set,  the  appropri- 
ateness  of  our  choice  is  evident.   It  is  also  immediately  clear  that  this 
set  of  vector  functions  actually  constitutes  a  linear  space  obeying  postu- 
lates (1),  (2),  and  (3)  if  to  the  abstract  +  operation  of  postulate  (l) 


■.;;.(  I'' ■ 


',  }    .     ■■<  ■■_■■ 


!'■.*.:»;.  i;  "I ; 


f    .'if.i   •  t~:  ■''■;i 
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we  make  correspond  the  operation  of  ordinary  vector  addition,  and  to  the 
abstract  nniltipllcatlon  operation  of  postulate  (2)  we  make  correspond  the 
multiplication  of  a  vector  "by  a  complex  scalar. 

To  coniplete  the  correspondence  between  the  abstract  space  and  its 
representation,  we  must  define  what  we  mean  by  inner  product,  linear  opera- 
tion, and  adjoint  operation  as  applied  to  the  vector  fxinctions.  Vedenote 
the  correspondence  between  any  abstract  entity  and  its  explicit  representa- 
tion by  an  arrow  -> .  Then  if  t.,->Z,j(z)  and  i^-^2,Az) ,  the  inner  product 
(t^.Cj)  is  defined  by  the  correspondence 


86) 


(Ci.C2)->/^^U*(:=)  •  y«)dz 


0 

where  the  asterisk  means  complex  conjugate  and  the  dot  indicatea  ordinary 
scalar  product  of  two  vectors.   It  is  trivial  to  verify  that  the  inner 
product  as  thus  defined  possesses  all  the  properties  attributed  to  the 
abstract  inner  product. 

Next  we  define  two  kinds  of  linear  operations  as  applied  to  the 
vector  ftmctions.  If  C->Z(z),  then  the  linear  operation  MC  is  defined  by 

87)  Mt  ->M(z)  Z(z) 

where  M(z)  is  a  square  matrix  of  rank  3  whose  elements  are  continuous  com- 
plex functions  of  z  over  the  interval  0  <^  z  ^  a.  Naturally  M(z)  Z(z)  indi- 
cates matrix  multiplication  with  Z(z)  regarded  as  a  column  matrix.  The 

linear  operation  N{,  on  the  other  hand,  we  define  &n    — 

a        _ 

88)  NC  ->y^  N(z,z')  Z(z')  clz« 

0 

\fhere  N(z,z*)  is  a  matrix  similar  to  M,  but  whose  elements  are  continuous 
complex  functions  of  both  z  and  z"  over  the  region  oSs^a.  0^z'Sa» 
except  for  a  possible  finite  Juflp  at  z  =  z*. 

It  is  equally  trivial  to  verify  that  the  operations  defined  by 
correspondences  (8?)  and  (88)  possess  the  defining  property  of  the  abstract 
linear  operation. 

We  define  the  operation  adjoint  to  M  by 

89)  M«C-»M«(z)  Z(z) 

where  the  matrix  M'(z)  is  obtained  from  M(z)  by  intercliange  of  rows  and 
columns  in  the  latter  and  replacement  of  all  elements  by  their  cosiplez 
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conjtigates.     We  shall  verify  that  if   t— ^V(z),  then 

90)  (M«:,t  )  »  (C.Mt)  -.      - 

in  accordance  with  the  defining  property  of  the  adjoint.  We  have 

91)  (M'C.f  )-»/  [m«(z)Z(z)]*  •?(z)dz  =y  Z  [ZM'ij(z)Zj(z)]*  V^(z)dz..^ 
»/^A  ^M«*(z)"z*(z)7.(z)dz  =/  2.  5_  Z*(z)M  (z)7  (z)dz 

-'o  i=i  j»i  ^    •'  /^     0  j=i  i=i  -^   •'^   ^     •:..;.:: 

=  /*5Izt(z)  r5lM,.(z)y.(z)  Idz  =/^Z*(z)  .  M(z)  V(z)  dz^(C.Mt) 
0  J=l  •'    •-  1=1  •'^    ^    •■    -^0 

which  proves  equation  (90). 

Finally,  the  operation  adjoint  to  N  is  defined  "by 

a 

92)  H'C  ->y^  N«(z.e«)Z(z)dz 

0 

where  the  matrix  N»(2iz')  hears  the  same  relationship  to  Il(z,z')  as  M'(z) 
does  to  M(z).  It  should  be  noted  that  the  integration  in  (92)  is  with 
respect  to  z.  Again  letting  t    ->•  Z(z)  and  £-»'7(z)  we  shall  verify  that 

93)  (N'Ca)  =  (CUE). 

&        fit 

94)  (N'5.t)  -->7'dz«yH<*(z,z«)  Z*(z)«V(z«)dz 


=  /'*dzt/*2I  21  H'*(z.z»)  Z*(z)V,(zOdz  =/dz»  /Xzt(z)  5lN..(z,z«)V.(z»)dz 
0    0  i=l  4=1  ^i  ^         ^  o     0  J=l  -^    i»l  ^^      ^ 

=r/'*Xz!(z)dz/*^II,,(z,z»)  Y.(z')dz«  =/dz  Z*(z)  7lT(2,z«)  ?(2')  dz»-XWt) 
^0  J=l  J     0  i=l  J^        ^         0  ^ 

where  the  reversal  of  the  order  of  integration  is  Justified  hecause  all  of  the_^ 
functions  are  continuous  except  for  the  possihle  finite  ^vaap  in  ^.^iz,z')eit  z=t^ 
Thus  equation  (93)  is  proved^ 


,'-fT  ::,Mvoi.^  V--'  '-^ '•'■';''■"?'■>'•  V''-u% 
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With  the  correspondence  between  the  abstract  linear  space  and 
the  linear  space  of  the  vector  functions  now  completed,  we  can  proceed 
directly  to  the  derivation  of  the  variational  fornrulas  for  E.  .  and  S. ,. 
To  1365:10  with,  consider  the  variational  form  hLC,C_]  girea  by  equation 
(73).  Let  us  permit  the  fixed  element  a  to  be  either  one  of  two  particular 
fixed  elements  a.  (i  =  1»2)  and  the  fixed  element  a*  to  be  either  one  of 

two  particular  fixed  elements  a.  (l  =  1,2).  Similarly,  let  the  variable 
element  C  "be  either  one  of  two  variable  elements  tj(j  =  1,2)  and  the 
variable  elemont  C  either  one  of  two  variable  elements  C.  (i  "  lj2). 
Introducing  these  possibilities  into  equation  (73)  we  get 

.     (L'C!,a  )(a'  LC.)  (j  «=  1,2) 

thereby  obtaining  four  different  variational  forms  of  the  type  defined  by 
equation  (73) »  corresponding  to  the  four  possible  subscript  pairs  i,j  (i  s  1,2), 
(j  =  1,2).  The  variational  property  of  hQ.  CrJ  may  evidently  be  expressed 

by  the  statement  that  hQ.,tr]  is  stationary  with  respect  to  small  variations 
of  Cj  and  Cj  about  the  elements')^,  and  ^.  respectively,  where  0^.  satisfies  the 

equation 

96)  aj  «  ZL^j  (J  =  1.2) 
and 01  •  the  equation 

97)  a[  »  K«  L"^^  (i  =  1,2) 

For  convenience,  let  us  denote  by  h.  .  the  value  of  ^[t^stj J  at  its 
stationary  point.  According  to  equation  (85),  we  have 

P      ,T    (L«/n«,a  )(a'  L/y  )     ,  ,     (j  =  1,2) 

98)  h,j  .  4^y^,  ]  '-j;;jP±^ ^=  K.  I'^j)  =  a'^i.aj)  (i  ,  1,2) 

We  now  make  two  observations.  The  first  is  that  linear  operations 
and  inner  products  of  the  type  defined  by  correspondences  (86), (87),  and  (88) 
appear  in  the  integral  equation  formulation  of  our  electromagnetic  problem 
and  in  the  integral  fonnulas  (66)  and  (72)  for  S. .  and  S.  ..  The  Becond, 
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which  we  shall  verify  below,  is  that  these  linear  operationg  together  with 
certain  electromagnetic  vectors  can  "be  put  into  correspondence  with  the 
abstract  operations  and  elements  occoxring  in  equation  (98)  in  such  fashion 
that  when  the  abstract  operations  and  elements  are  replaced  by  their  electro- 
magnetic counterparts,  the  resulting  value  of  h. .  is  H. .  or  S. .,  dependii^ 
on  the  precise  correspondence  chosen. 

Ve  vei'ify  the  above  assertion  first  for  E  . ,  obtaining  for  the 

latter  the  desired  variational  formula  as  the  end  result  of  the  verifica-  

tion.  The  necessary  correspondences  are  set  up  in  a  systematic  way  as 
follows:  starting  with  the  integral  formula  (66)  for  H^^.,  which  we  repro- 
duce here 

0 

we  recognize  that  according  to  the  definition  of  inner  product  given  by  cor- 
respondence (86) ,  we  shall  have 

100)  (a^  ,  LOjj)  =  E^j     ■  ' 
if  the  following  correspondences  are  made; 

101)  a[   ^u;e-^^«°««^ 

103)  v^-^hi)^^^ 

How  from  its  definition,  E/.v(z)  satisfies  the  integral  equation 


10/.)  S.  e"^'   ^»««-  [l^cM]\^^U)'J^f-^''^''''^''''^^^*-)^i-%f-^ 


(J  =  1.2) 
Therefore,  in  view  of  correspondence  (103).  we  can  mak»  equations 
(96)  correspond  to  equation  (104)  if 

,--»  ^-  ^ikz  cos* 

105)         .  a,   — >u.  e 
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and  if 


0 

whence,   from  correspondence   (102)  we  olJtain 

a 

107)  E->Cl+C(z):    2_S^B-l(^j_y-^ik|.-..|    C08<)^(^^^,j    ^^^,^^, 

There  remains  to  establish  a  correspondence  'between'^,  and  some 
yector.   Since ''^.  satisfies  equation  (97),   we  set  up  the  integral  equation 
corresponding  to  equation  (97)  and  define  the  vector  corresponding  to^. 
as  the  solution  of  this  integral  equation. 

Tor  this  purpose,  we  must  form  the  adjoint  operations  L*  and  K' 
according  to  the  defining  correspondences  (89)  and  (92). 

From  correspondence  (102)  we  have 

108)  L.  ^  -  j^  B.(z) 

while  correspondence  (107)  gives  -        " 

109)  K'  ->  -^^Br^ri*C.(z)D  -/V^I^-^'l°°««A(z.,0...dz»  * 

where  to  obtain  the  first  term  on  the  right  of  (109)  from  the  corresponding 
term  in  (107)  we  have  used  the  fact  that  (B"-*-)  =  (B*)   for  any  non-oingular 
matrix  B  plus  the  fact  that  if  7   and  G  are  any  two  matrices,  than  (?G)  -^  G  T  ti — -  ■ 
Combining  correspondences  (108)  and  (109),  w«  get 

110)  z'l'->[i  +  B'-^(z)  C.(z)B«(z)  ]  +  j^/V"^l«-^'I°°^*A(z«,z)B'(z0...dz« 

o 
The  first  term  on  the  right  of  (llO)  can  be  simplified.  Trom  the 

definition  of  c(z),  we  see  that  C'(z)  is  a  matrix  whose  first  two  columns  are 

zero  and  whose  last  column  equals  that  of  B'(z)«  Therefore,  we  have  at  once 

0  0  o' 
(111)  B«"^(z)C«(z)   »   I  0  0  0 

0  0  li 


*  Ai(s'*z)  "  A(z*,s)  since  tha  matrix  A  is  real  and  symmetric. 
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Trom  eqtiatlon  (111)  there  follows 


.  /000\  /     0       0     0 

112)  ^*      M  0«(z)  B»(r)  »     000     B«(z)-        Q       0     0 

\,00l/  ^B'3,    B«2  B^3 

/    0       0       0 
and  since  C(z)  was  defined  as  the  matrix         0       0       0      |  ,   It  is  clear 

that  the  last  term  on  the  right  of  equation  (112)  "bears  the  same  relation- 
ship to  B'(z)  as  0(z)  does  to  B(z).     Denoting  this  last  term  by  0(z),   the 
final  form  of  K  L    "becomes 

113)  ZL     -^    [l  +  ^(z)]+     ^^yVi'^>^-^'i°°««A(z..z)B.(z.)    ....dzt 


Now  making  iise  of  correspondences  (97),  (101),  and  (113),  and  denot- 
ing "by  ^^i^)   the  vector  obeying  the  correspondence 


_» 


m)  ^i  ->l^j(z) 

it  follows  that  ^A^)  satisfies  the  integral  eq.iaatlon 

0 

(i  =  1.2) 

which  has  a  form  precisely  analogous  to  that  of  the  integral  equation  (lO'l). 
Equation  (115)  is  called  the  %dJolnt«  equation  to  (104).   Its  physical  inter- 
pretation will  "be  given  in  Section  6,  which  deals  with  the  reciprocity  theorem 
referred  to  in  the  Introduction. 

We  have  thus  set  up  the  electromagnetic  counterparts  of  each  of  the 
abstract  elements  and  operations  appearing  in  equation  (93),  as  well  as  the 
counterparts  of  equations  (96)  and  (97).  Combining  equations  (98)  and  (lOO) 
we  have  ,.-•?.'■'■'. 

116)  E..  -  7 '^ -^ 

'       ^J      (^.,  LKL'jj,)     :- -— 
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and  replacing  each  of  the  ahstract  entities  In  formula  (116)  by  Its  electro- 
magnetic counterpart,  we  obtain  Immediately  the  final  yarlatlonal  formala 
for  R, 


ij 
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Because  of  the  stationary  character  of  the  form  ii^t-.tij  and. 
the  correspondence  we  have  established  between  the  abstract  space  and  the 
space  of  vector  ftmctions,  we  can  siun  up  the  properties  of  formula  (117) 
in  the  following  way:   if  Z.  x(z)  satisfies  the  integral  equation  (104) 

and  E/  v(z)  satisfies  the  adjoint  equation  (115),  then  of  course  formula   . 

(117)  gives  the  exact  value  of  R.  .;  on  the  other  hand,  if  for  E/.x(z)  and 
E/  x(z)  we  substitute  approximate  trial  functions  into  formula  (117),  then 
the  formula  will  be  stationary  with  respect  to  snaall  departures  of  the  trial 
functions  from  the  true  solutions  to  the  integral  equations  (10^)  and  (115), 
respectively.  Naturally,  it  is  the  second  property  which  makes  formula  (117) 
more  accurate  than  formula  (99)  for  the  purpose  of  approximate  calculation 

of  E^j. 

The  procedure  for  derivir^  the  variational  formulas  for  S^,  is 
exactly  like  the  preceding  one  save  that  the  electromagnetic  counterpart* 
of  a  and  7).  mast  be  changed.  Eeproducing  here  the  integral  formula  (72) 

0 

we  see  that  we  shall  have 

119)  (a'.  I'^j)'  S^j  ■ 

if 

•  ^_   ikz  cobQ 

120)  ttj^  "^^i  ® 

while  correspondences  (102)  and  (103)  remain  unchanged. 

If  we  again  define  E/  v(z)  by  the  correspondence 

la)         ^;-»e[^^(z) 

where  ^  '  is  the  solution  to  equation  (97) .  then  it  follows  from  equation 
(97)  and  correspondences  (113)  and  (120)  that  Ej^(z)  now  satisfies  the  inte- 
gral equation 

Eqtuatlon  (122)  is  the  adjoint  to  eq.uation  (104).  Th\is  the  adjoint 
integral  equation  associated  with  S^^  differs  from  that  associated  with  B^j, 
viz.,  equation  (115).  only  in  the  "incident  wave"  term  on  the  left.  The  phys- 
ical interpretation  of  equation  (122)  will  also  be  given  in  Section  6. 
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6.  The  Reciprocity  Theorem 

In  the  prerioua  Bection  we  set  up  certain  correspondences  between 
the  abstract  entities  of  the  linear  space  and  the  vectors  and  operations  of 
the  electromagnetic  problem,  and  we  showed  that  the  correspondences  could--- 
be  80  chosen  that  the  inner  product  (a^,  L  7]^)  equaled  E^^  °^  ^^y   wiiicherer 
was  desired.  We  did  not,  however,  inquire  into  the  resulting  electromag- 
netic significance  of  the  analogous  inner  product  (L  1^^,  a^   and  we  wish 
to  do  80  now.  Considering  first  the  case  where  the  correspondences  are 
those  associated  with  R^j,  we  shall  show  that  (L  ^  ^fO.^)   can  be  inter- 
preted as  the  Ji*^  element  of  the  reflection  matrix  in  a  new  electromag- 
netic problem  whose  relationship  to  the  original  problem  is  a  very  eiiaple 
one.  We  ahall  show  farther  that  one  of  the  equalities  in  equation  (98), 
namely, 

126)  (aj  .  LTl^)  '   (L*7?i.  aj) 

is  actually  the  expression  of  an  interesting  reciprocity  theorem  which 
connects  the  reflection  matrix  of  the  new  or  •reciprocal"  problem  with  the 
reflection  matrix  of-  the  original  problem. 

We  deduce  the  statement  of  the  reciprocal  problem  in  the  follow- 
ing way:  according  to  correspondences  (105),  (108),  and  (114),  together 
with  definition  (86)  of  the  inner  product,  we  have 

127)  (L  V^   .  aj)  -  27^  /  **        ^J^  ^^   (i)^^ 
where  b!  At)   satisfies  the  adjoint  Integral  equation  (115)^ 

A  comparison  of  expression  (99)  with  the  right  side  of  equation 
(127)  suggests  that  the  latter  may  be  interpreted  es  the  Jl   element  of 
the"  reflection  matrix  for  a  reciprocal  problem  In  which  the  scattering 
medium  has  a  susceptibility  matrix  b'*(z)  while  the  electric  field  within 
the  medium  is  l!*j(z).  The  exact  statement  of  the ^reciprocal  problem  is 
inferred  from  the  integral  equation  satisfied  by  B(^j(z).  Since  S(^^(z) 
is  the  solution  to  the  adjoint  integral  equation  (115).  It  follows  from 


the  linearity  of  this  equation  that  E^j^(b)  satisfies  the  complex  conjugate 


U-"^  ■■■■■■.::■  :H 
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of  the  adjoint  Integra.!  equation.  I.e., 

(1  =  1.2) 

This  le  the  integral  equation  for  the  reciprocal  problem.  To  inter- 
pret it,  we  compare  it  with  the  integral  equation  for  the  original  problem, 

TlZ., 

0 

(J  =  1.2) 

The  first  point  of  difference  to  be  noted  between  equations  (128) 
and  (129)  is  in  the  susceptibility  matrices.   It  will  be  recalled  that  B»(2) 
was  derived  from  8(2)  by  transposition  of  rows  and  columns  in  the  latter  and 
replacement  of  all  elements  by  their  complex  conjugates;  B**(2)  is  therefore 
just  the  ordinary  transpose  of  B(z),  Moreover  C(z)  was  defined  as  having  its 
first  two  rows  zero  and  its  third  row  equal  to  the  conjugate  of  the  third  column 
of  B(z);   therefore  %   (z)  bears  the  identical  relationship  to  B**(z)  as  C(z) 
does  to  B(z).  We  conclude  that  in  the  reciprocal  problem,  the  scattering 
medium  has  a  susceptibility  matrix  which  is  the  transpose  of  the  matrix 
characterizing  the  original  medium.  We  call  the  new  medium  the  "transposed 
medium  ■. 

The  second  point  of  difference  between  the  two  Integral  equations 
is  that  the  matrix  A(z,r«)  in  equation  (129)  has  been  replaced  by  A(z«,z)  in 
3quation  (128),  Eeferring  to  equations  (26)  and  (27),  which  define  A(z,z«) 
fe   see  that  interchanging  z  and  z»  is  equivalent  to  stibstituting  -€  for  4 
In   the  matrix.   Inspection  of  the  remaining  Q  terms  in  equations  (128)  and  (129) 
.;hen  shows  that  as  far  as  the  «  dependence  is  concerned,  equation  (120)  is 
)btainable  from  (129)  simply  by  replacing  «  in  the  latter  by  -«  ,  We  con- 
Jlude  that  in  the  reciprocal  problem,  the  angle  of  incidence  is  -«  . 


".  i^  .f  .^'.  'Z  K  -4  v*v  • '' 
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Finally,  it  ifl  clear  from  eqiaation  (128)  that  the  Incident  wave 
vector  amplitude  for  the  reciprocal  protlem  is  u.  .  This  observation  com- 
pletes the  definition  of  the  reciprocal  problem. 

Nov;  in  Section  U,   we  derived  the  integral  fonnulfr  (66)  or  (99) 

for  E.  .  The  same  derivation  obviously  applies  at  once  to  the  reciprocal 
xj 

problem;  in  fact  if  we  denote  by  E'  the  reflection  matrix  for  the  recipro- 
cal problem  we  have 


0 

Combining  equations  (100),  (126),  (127),  and  (130),  there  results 
131)  E,j  =  e'^ 

Equation  (131)  constitutes  the  reciprocity  theorem  for  the  reflec- 
tion matrix.  The  verbal  statement  of  this  theorem  is  68  follows:   in  the 
original  problem,  let  an  incident  wave  of  vector  amplitude  u.  and  angle  of 
incidence  ft  excite  a  reflected  wave  whose  amplitude  component  in  the  u. 

direction  is  R. .;  then  an  incident  wave  of  vector  amplitude  u^  which  im- 
pinges  on  the  transposed  medium  at  an  ajigle  of  incidence  -0  ,  will  excite 
a  reflected  wave  whose  amplitude  component  in  the  u.  direction  is  also  equal 
toR.j. 

The  reciprocity  theorem  for  the  transmission  matrix  is  derived  in 
similar  fashion.  Using  correspondences  (105),  (108),  and  (121),  we  obtain 

"3)    (L--,  ;  .  o,)  =  ^,  />^   -=*u^-  b"(0  f-  CO  d. 


0 


where  E,  v(z)  now  satisfies  the  adjoint  equation  (122).  From  the  form  of 
expression  (132),  it  is  again  reasonable  to  define  the  reciprocal  problem 
as  being  governed  by  the  complex  conjugate  of  the  adjoint  equation.  Thus    v 
the  integral  equation  for  the  reciprocal  problem  is  now 

o 

For  comparison,  we  also  write  the  equation  of  the  original  problem 

once  more 

134)  u/'''"'*-[l*C(Oli(j)(0  -  2~rc  /e^l^-^'l*'"^^A(z,z.)B(«.)l^j)(z.)dz. 
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From  eqiaation  (133)  it  Is  clear  that  In  the  transuiBsion  matrix   "^ 
reciprocal  problem  the  original  Bcatterijag  medium  la  again  to  "be   replaced 
"by  the  transposed  medliam.  However,  the  complete  epeclflcAtlon  of  the  re- 
ciprocal problem  requires  slightly  more  detailed  consideration  of  the  con- 
Jugate  adjoint  equation  than  was  necessary  In  the  case  of  the  reflection 
matrix. 

Conalder  first  the  form  of  E^  »(z)  in  the  region  s  >  a.  Since 
in  this  region^  (z)  ■  0  and  |z-2*|  =  z-z',  equation  (133)  gires 

T»c\  5**  /   \       -  ^-ilcz  cos*,  ^ikz  COB 4   Ik  r-ikz'coB©  .  _•♦,  ^•♦/  .\.    . 
135)  B^^j(z)»u^e        +6       i^^To/*        V  (^^(i)(^')4«» 

z  >  a 

-  I*  .  . 

Thus  in  the  region  z  >  a,  ^(^\M   represents  a  downward  moving 

Incident  wave  of  amplitude  u.  impinging  on  the  scattering  medium  from  above, 
plus  an  upward  moving  reflected  wave.  From  the  fact  that  A(z,z')  in  equa- 
tion (13*1)  is  replaced  by  A(z«,z)  in  equation  (133),  it  is  easily  inferred 
that  the  x  dependence  of  the  field  vectors  In  the  reciprocal  problem  mast 
be  e"   3C  aln«  ^  ^^  conclude  that  in  the  transmission  matrix  reciprocal 
problem,  the  Incident  wave  strikes  the  transposed  medium  from  above  at  an 
angle  of  incidence  tt  +  0-  with  the  positive  z  direction,  where  Q   is  the  angle 
of  incidence  in  the  original  problem. 

In  the  region  z  <  0,  C  (z)  =  0  and  |«-z«l  «  ««-b.  Equation  (133) 
then  becomes 

ft 


=  »*    /   X       -     -ikz  co8«-      -ikz  cos  Q    Ik         /rlkz'cosC.  ,,*•/,  vs'*  (,,> 
E^^)(2)  -  u^e  +e  5^^^  ye  A_^B  (z')E^j(z«j 

o 

z  <  0 

.1* 


Therefore  in  the  region  z  <  0,  E/.x(z)  represents  a  downward  moving 
transmitted  wave  whose  propagation  direction  makes  an  angle  v  +  Q  with  the 
positive  z  direction. 

Now  let  T*  denote  the  transmission  matrix  for  the  reciprocal  prob- 
lem. Taking  the  scalar  product  of  equation  (136)  with  u.  and  dropping  the 
factor  «"*    °°®  ,  It  follows  from  the  definition  of  tranamlBsioB  matrix 
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eirea   in  Section  (h)   that 

which  in  view  of  equation  (63)  and  the  BTminetry  of  A  ,  reduces  to 

If  we  define  the  matrix  S'  throtigh  the  relation 

139)  t'  "  I  +  SI 

we  have  from  equation  (138) 

Ikz'^cosA-    _'*/.\='*  /  .\  ,  . 
e  *    ^j  '  B  («')  ^(i)(z^  ^* 


lilO) 

^Ji  = 

ik 

2C0B« 

I 

0 

ComlJinlng  equations 

141) 

hr 

or  equivalentl7 

142) 

'ir 

t' 

Bqoatlon  (1^2)  espresses  the  reciprocity  theorem  for  the  trans- 

mieelon  matrix.   In  words  it  may  "be  stated  as  follows:  in  the  original 

_  and 
problem,  let  an  incident  wave  of  vector  amplitude  u./ angle  of  incidence  0 

J 

excite  a  trensmitted  wave  whose  amplitude  component  in  the  u.  direction  is 
T. jj   then  an  incident  wave  of  vector  amplitude  u.  impinging  from  above  on 
the  transposed  medium  at  an  angle  of  incidence  tt  +  9  with  thfi  positive  z 
direction,  will  excite  a  transmitted  wave  in  the  region  z  <  0  whose  ampli- 
tude component  in  the  u.  direction  is  also  equal  to  T.  .. 

The  relationship  between  the  original  problem  and  the  two  recip- 
rocal problems  is  indicated  graphically  in  rig«  2.  The  lines  with  arrows^ 
represent  the  propagation  directions  of  the  resi>ective  waves.  ,   r.  r~~ 

+  :\-v.:!:rf:;';n  "' '  ', ,    _,.:;*.•:';"':':•;¥;;■  vir^^^;  .■■;;  -.j.^i^^-;  ■■■^..  -  •:  •- .M' ,.,jV'.'X-'<  ^-.;^^  .'""-v^  ,  j- 

'■^Trr,-.:.!:?  ^^-^  .;t"J>  .?'i?■i:^'^:t■.  ■'  lOc*  ^^n- ,'.*  ■^'•:>■■■'^  t  ?■  Tie  ?  ;.f-"~ '•";■  '^■'"f    ■;-'^"?   .  n;.  .■•♦i:^  T" 
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Reflected  Wave      Incident  vrave 
Reciprocal  Problem  for  Reflection  Matrix 


Reflected  Wave 


Transmitted  Wave 

Reciprocal  Problem  for  TranamiflBion  Matrix 


rig.  2 

We  conclude  with  an  example  which  will  illustrate  the  physical 
relation  between  the  original  scattering  medium  and  the  transposed  medium. 
Assume  that  the  region  0  5  z  ^  a  is  filled  with  a  mixture  of  electrons  and 
positive  ions.  The  electron  and  ion  densities  are  equal  at  every  point, 
depend  only  on  z,  and  naturally  vanish  at  z  =  0  and  z  =  a.  We  assume  that 
there  is  also  present  in  the  region  some  magnetic  field  H  which  is  in 
general  a  function  of  z  but  constant  in  time.  Finally,  we  suppose  that 
the  mass  of  an  ion  is  very  much  greater  than  that  of  an  electron.  Such 
a  mediim  (known  as  a  magneto-icnic  meditun)  is  generally  used  to  represent 
the  earth's  ionosphere,  H  then  being  the  earth's  magnetic  field. 

Now  let  IT(z)  be  the  electron  or  ion  density  function,  e  and  m 
respectively  the  charge  and  mass  of  the  electron,  ui  the  angular  frequency 
of  a  time  harmonic  electromagnetic  wave  impinging  on  the  medium,  and  »(z) 

the  number  of  collisions  per  second  which  a  given  electron  makes  with  Ions 

22 

or  atoms  .  Then  the  susceptibility  matrix  of  such  a  medium  can  be  shown 

22.  Wilkes, M. 7.  Proc.Roy.Soc.  A,  vol.  189,  p.  130,  19^7.  See  equation  (2.2) 
of  Wilkes'  paper. 
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to  be  /  \ 

Y^  -d+iZ)^         iY^(l+iZ)+Y^y2     Yj^Y^-iYgd+iZ)  ' 


1^3) 

B(0  ^ 


(l+lZ)[(l+iZ)24r2+Y2+Y^ 


YgTj^  -lY^(l+iZ)         y|  -   (1  +1Z)^     lYj^d+iZHYgY^ 
^1Y2(1+1Z)+Y^Y3^  Y^Y^-iY^ftHZ)       Y^  -(l+iZ)^ 


wr,^«2  eHe  (z) 

where  X(0  =  \   ^  .       T.  (O  =     T^ (i  =  1.2.3) 

0 

and  Z(8)  =     ^^^'      .  H  is   the  x,y,   or  z  component  of  H 


^     .^^  .„  .^.  -..,.  ..   ^ -g 


Inspection  of  equation  (1^3)  shows  at  once  that  transposition  of  the 
matrix  B(z)  is  equivalent  to  replacing  Y^^,  Y^,   and  Y-j  hy  -Y^^,  -Y^  and  -Y^ 
respectively.  Therefore  if  the  original  scattering  medium  is  a  magneto- 
ionic  medium  of  the  type  Just  defined,  the  transposed  medium  differs  from 
the  original  only  in  that  the  magnetic  field  H^  Is  reversed  in  direction 
at  every  point. 


7.  Conclusion 

Schwinger's  integral  equation-variational  technique  has  been 
applied  to  the  problem  of  the  scattering  of  an  elliptically  polarized 
plane  electromagnetic  wave  by  a  continuously  stratified  anisotropic  medium. 
Variational  formulas  have  been  obtained  for  the  elements  of  the  reflection 
and  transmission  matrices  of  the  problem.  The  efficacy  of  certain  consider- 
ations in  the  theory  of  linear  spaces  with  respect  to  construction  of  the 
variational  expressions  has  been  demonstrated.  ^Finally,  a  certain  symmetry 
property  of  the  variational  formulas  has  been  used  to  prove  a  reciprocity 
theorem.     -fi  v---  '  -    ■■■<■-■■  ■'.:-:\'-:h.   '■  ..•''^9>.lvr^''- '  ;^''i.:;^;^iv  v'J.j;-?  :!.  ;>;  ■'■■-   .  .^  ,  . 


'..yAfS  ^.v-^^.^^-.  ,-^h-:\    U'r-^)  'b^'f   T-:'-')  "^s :i.? :^^;;•0  Ci 7' 


I,       '-•   c . 


i'l 


jftr  .'do: 
t  ■  ♦■  r 


:..\  .. 
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Appendix  I 

We  prove  that  If  f  and  g  are  solutions  to  Mazvrell's  equations 
in  a  region  where  M-  «  p.^  and  e  is  a  continuotis  "but  otherwise  arhitrar7 
matrix  function  of  position,  i.e., 
la)  V*f  -  Icup-^g  =  0  ^^-—^ 

2a)  VXg+lu.cf-0, 

then  there  always  exists  at  least  one  rector  function  ^  satisfying  the 
equation 

»  V^^  +  k^^  =  »  ^  (e  -  c^X)f 

o 

in  the  given  region,  and  having  the  further  property  that 

^)  f  »  7(7  •^)  +  k^^ 

5a)  i  =  -  i(jue  7x^ 

0    ' 

Taking  the  divergence  of  eq\aatlon  (1«),  we  get 
6a)  7  .  g  a  0 

Since  g  is  continuous  with  continuous  first  derivatives  in  the 

region  where  it  satisfies  equations  (la)  and  (2a),  then  l)y  a  fundamental 

theorem  in  vector  analysis  there  exists  a  vector  function  i>      such  that 

•  0 

7a)  g  B  _ioj«  7y^ 

0       0 

in  the  given  region,  where  the  constant  factor  -iwe  is  added  for  con- 
venience. 

Substitution  frdm  equation  (7a)  into  equation  (la)  gives 

8a)  7  X  (f  -  k^<6  )  =  0 

-    2  - 

Thus  the  vector  f  -  k  ^  ^  is  irrotational  in  the  given  region. 

By  another  fundamental  theorem  in  vector  analysis  there  mast  exist  a  scalar 
function  q^  such  that 

9a)    f  -  k^<5   =  -  7q 

J.,  Let  us  write  equation  (2a)  as 

10a)  7Xg  +  ia,c  f  =  -iu»  (c  -  e  I)  f 

0  0 

Suhstltuting  from  equations   (7a)  and  (9a)  into  equation  (10a),  ve 
obtain 

11a)  -ia.<^7x7)t^^  -    iuic^  ^^+4u«^k^<f  ^  •    -l(ii(c  -  c   l)f 
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which  after  dlylsion  "by   i«juc    becomes  

0  _   -  .  -      - 

12a)  -7XVX^    -   TTqt    +  k^A     =  -     -i  (e  -  c   I)  f 

•  0  0  •  0  C  0 

o 
Equation   (12b)  is  a  necessary  condition  on  any  pair  of  functions 

<p     and  Q     which  are  related  to  f  and  g  by  eqiiationa    (7a)  and  (9a), 

Now  let  <p    and  a     represent  one  partictdar  pair  of  frinctiona 
satisfying  equations   (7a),    (9a),  and  (12b);     then  the  functions  <^  and  q 
defined  "by 

13a)  ^  =  ^jj  +  Vp 

1^)  q     o     q     ■«•  k     p 

where  p  is  an  arbitrary  diff erentiable  scalar,  also   satisfy  equations   (7a), 
(9a),  and  therefore   (12b),   i.e., 

15a)  i  B     -iujc     7X^  • 

o 

I6a)  f  »    -  7q  +  k^^ 

17a)  -  7X  7xf  -  7q  +  k^<f «  -     ^  (e  -  e^l)  f 

0 

Since  p  is  arbitrary,   let  us  choose  it  to  be  a  solution  of  the 
eqtiation 

18a)  7^  p  +  k^  p  =  -7  '^^  -  q^  ' 

•  o        o 

Equation  (I8a)  may  be  wtitten  as  •■•..( 

19a)  V  .  (^^  +  7p)»  -  (q^  +  k^p) 

or  .   ' 

20a)  7  '^s  -  q 

in  view  of  equations  (13b)  and  (1^).  ' 

Substitution  of  the  value  of  q  given  by  equation  (20a)  into 
equation  (17a)  gives  "'  - 

.2J:  _    1. 

c 

o 

which  in  rectangular  coordinates  becomes  ,       1  ■  ^ 


21a)  -.7X7)C^  •♦■  7(7  '^  )  +  k^^  =  -     =-  (c  -  €^l)  f 


/  '■. 


22b)  T^f  +  k^^  =  -     7"  (c  -  €qI)  f 

"-       Under  the  same  substitution,  equation  (l6a)  becomes 
23a)  _^ ^ f  =  7(7  '(f  )  +  k^<f     .  ;     ' 


-  k2  - 


Equfttione  (22a),  (23a).  and  (15a)  agree  with  equatione  (3a), (J^), 
and  (5a)  respectively.  Therefore  the  function  4*  defined  "by  equation  (13a), 
where  p  satisfies  equation  (18a),  oheys  equations  (3a),  (^),  and  (5a).  The 
existence  of  the  desired  (p   is  thereby  proved. 

Appendix  II 

In  this  appendix    carry  out  the  evaluation  of  the  term 

ikxsLnO 

k  cos '_  .^ 
o       "^0 

(see  equation  2^),  where  the  differential  operator  in  ea^aression  (It)  appll 
of  course  to  the  variables  x  and  z.     Beginning  with  the  operation 

(7  •    )  "     i     ^~  +  k    r— and  ignoring  the  constant  term    -rr — TTTTrZI^®  have 

~*      OZ  ""      OZ  fcX€    K    COB  ^ 


Ih)    t(V  .    )  +  k!l     e^^^"^^        /  ^lk|z-2«lcos  «.  j(     J 
'    ^  ^  -•    2i€  k  coefi     /  ^  • 


es 


0 


2b)    (i  I-+  k  I-)    •  s^^  8inOy-\ik|z-zMcosfl^(^,j  ^^, 

—  dx       —  dz  -f 

0 

.   (i  |_+  k  |_).ei^^  «i««  r/e"^(^-^'^'^°"%(z')dz'  +ye^^(^'-'^°°«%(z')dz' 1 
The  operations  on  the  right  of  equation   (2b)  result  in  ,, 

3b)      i     1^.    p^«^«_/V'^(^-^')*^«<»P(z«)dzt] 

=  ik  8in;>e^^  Bin«y^lk(z-z')co8  Q  p^^,,)^,, 

0 

i,T,)     1    L..   Felfcr  8in«y*^lk(z«-z)cos«j(^,j^  n  .  :  : 

.,      .    -    ikx  sinfi  r  ik(z'-s)coB^  ^  /_,x   ._| 
«  ik  sinO  e  V  *  PiC^*)  *»' 

2 

5t)     k|_  .   Life.  Bln«y^ik(z-zt)cos«p(^,j  ^^,1 

o 


•  P,  «nd  P-  are  the  x  and  z  cooponenta  of  P,  respectively. 


-43- 


6b) 


k    |_.   pkx9ln«y^^ik(z«-.)co80  3^(^,j^^,-j 


.,  rt    ikx  BlnO/^  ik(z*-z)coB«  _  ^    ,v.    ,  ^Ikx  sin^  „  /   \ 

■  -  ik  coaQe  2J  ^  P~(s')dz'-e  PgU) 

z        • 
AddiDf  these  equations,  we  get 

7b)     7  -p^  «i^«yei^l^-='l°°««P(z.)dz'] 


s  ik  einO  e 


r^^lk(z-z')cOB;^^(^,j^^,^^ik(z'-z)c08«p_^(^,^^ 

*  ik  cosOe^^  ^^^«  [7'V^(^-='^°°^«P3(z')dz./e^^^^'-^)"«^z.Hz^ 

Next,  we  inust  take  the  gradient  of  the  right  side  of  equation   (7b), 
i.e.,  we  must  operate  on  it  with   (l.  |^  +  k  |^).     Tlxe  x  differentiation  gives 


after  a  little  rearranging 


'■      8b)     kV^  8in«y-%ik(z-z')co8«  T^^^a^  p^^^.j  _  smo  cob«  P^Cz')]    i  dz" 

BinY%ik(z'-«)cose  [l8i„2^p^(2,)+  BinOco9fiP3(z')]    1  dz' 


+  k2  e"^ 

z 

while  from  the  z  differentiation  there  results 


9b)  k^  e^^^  8in«y%il^(^=')<^o««  [Isind  eoeO  P^(z«)-co82«  P^Cz')]  ^4^, 

ik  e^^  ^  ®^^  r8inOP^(z)k+  coe^P^Cz)  kj 
+kV^  siny^ik(z'-z)co8«  LaOcosOPiCzO-cos^O  ^^i^')\  k  dzt 

+  ik  e^  *^^  IlsinOP^Cz)  k  +  co9«P3(z)  ^      "^^  ' -^ 


o 

+ 


z 


-  w»  - 


The  complete  operation  U(V  •    )  -♦■  k^       is  now  formed  1)7  adding  to 
the  sum  of  expreeeions   (ffb)  and(91>)   the  term 


101>)     ic2e^'^^«''^V«"'''""'"'''^(^')^^* 

o 

=  j^2  ^ikx  8ln«y^ik(«-8«)co8«  fp^Czt)!  +  PgCe'U*  P3(2')  k]   iz* 
o 

g;    7  +  k^  e^^  8inOy^lk(7«-Oco8«  Tp^c^i)!  +  P^Cz'U*  PjC^O  kj  d 

The  result  of  this  addition  may  "be  arranged  as  follows: 
m)     kV^  Bin«y;ik(z-s«)co8«  r   ^2^  p^(^,)  ^  o.P2(z«)-8in«co8fiP3(z')j   dz«   i 

+  k^  e^^  X8inOy^ik(z-s')cos«  ["o.p^Czi)  +  P^Cz')  +  O-P^Cz')]  dz«  J. 
o 

^  ^2  ^ikx  sin0y-%ik(z-z«)co8<J  f  gm^  cosO  P^(z«)  +  0'-p^iz*)^BiiihT^iz')jdz*  k 

0 

2     ltac8in«/lk(z»-z)co84"co82«P  (z«)+  O.P  (z')  +  BinOcoseP3(z')]dz«  i 
+  ke  _/e  L  -^ 

^  j^2  ^ikx8ln«y^ik(z«-z)  cosO  fo.p^Cr*)  +  PgC*')  +  O-P^Cz')]  dz'     A 

z  '- 

^  ^2  ^ikxBLnOy^lk(z«-z)  cobO  Tm;,  ,o8  0  Pj^Cz')  +  Q'-B^i^^y^^in^^^^U')^^^'  k 


+  2ik  9^^^"*co8«P3(z)  k 


,2    tkxsinO  /  lk|z-s«|co8«A(z.««)  P(z«)  dz'  *  2ik  e"=  *  «^^*  cob«P3(z)  k 
=  ke  ._/o  


where  the  mtrlx  A(z.z«)  Is  defined  ^7  e<iuationfl  (26)  and  (27). 


I 
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Plnally  we  maltiply  equation  (ll"b)  "by     oic  k  coa^      ***  oTjtaln 


*o  °®  O  0 

which  ia  eqviation  (25). 


■Vv      ,-->=*S  ;. 


*  ;> 


'■"•'••■    tie 


-  if6  - 
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